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Abstract 



We study the 3d Ising spin glass with ±1 couplings. We introduce 
a modified local action. We use finite size scaling techniques and very 
large lattice simulations. We find that our data are compatible both 
with a finite T transition and with a T = singularity of an unusual 
type. 
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1 Introduction 



Three dimensional spin glasses [1] are a fascinating subject. Numerical simu- 
lations are here particularly interesting [2, 3, 4, 5, 6, 7], since for such model 
(the real thing) it is very difficult to obtain reliable analytical results. Up to 
date numerical simulations for the Ising case have shown a phenomenology 
very similar to the experiments on real spin glasses (for recent simulations 
and analytical results about, for example, aging phenomena, see [8]). The 
study of small size systems (up to a linear size L = 14) has shown a reason- 
able agreement with the predictions of broken replica theory, but it is obscure 
how much information about the thermodynamic limit can be inferred from 
the behavior of small systems. In particular one has to be careful about 
extrapolating the pattern of replica symmetry breaking from small to large 
lattices. Here our aim has been to reconsider the whole subject and try to 
clarify the emerging physical picture at low temperature T. 

We will deal with the problem of the nature and the existence of a phase 
transition. A cursory look at the history of the subject is useful. If we 
look at the period that begins when people investigated first the subject of 
disordered spin systems we can easily establish that there have been periodic 
oscillations, with periods of the order of 7 years. Researchers in the field 
have been oscillating between the credence that there is a sharp transition 
and the belief that there are no transitions at all (as is maybe true in real 
glasses) and that when lowering T there is only a gradual freezing of the 
dynamical degrees of freedom. At the beginning theoreticians had (at equal 
time) a different credence from the experimental researchers. The two groups 
had a different frequency of oscillations, and now there is consensus that the 
system undergoes some kind of phase transition. 

We have run long numerical simulations at various temperatures, doing 
our best to distinguish among these two possibilities. In our analysis we 
have been very much inspired by the approach and the doubts of Bhatt, 
Morgenstern, Ogielsky and Young [2, 3, 4, 5, 6, 7], and over their results we 
have tried to build and improve. We have found that the whole set of our 
data is well compatible with the possibility that there is a transition at a 
given non-zero T. Such a transition would be characterized by a large value 
of the exponent 7, close to 2.5 (7 is the usual susceptibility exponent, which 
will be defined later in a more precise way). The whole set of data is also 
compatible with a large set of possible reasonable functional dependencies, 
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which imply a transition temperature of zero. Recent studies using improved 
Monte Carlo techniques [9, 10, 11] also find that doubts about the existence 
of a finite T critical behavior are justified [12]. The difficulty of resolving 
among the two behaviors is because a large value of 7 implies that the system 
is not far from being at its lower critical dimension (at which, according to 
the conventional wisdom, 7 — > 00). The distinction among a system at the 
critical dimension and a system very close to it is particularly difficult to pin. 
We believe however that we are not too far from being able to resolve among 
the two models and that an increase in the simulation time of one or two 
orders of magnitude can clarify the situation. Of course precise theoretical 
predictions on the behavior of spin glasses at the lower critical dimensions 
would be invaluable. 

We have been studying the 3d Ising spin glass, with ±1 couplings, but 
we have not used the standard first neighbor model. Hoping for some gain 
we have simulated a slightly modified model with second nearest and third 
nearest coupling. The reason for introducing this model is that in the con- 
ventional model (on the usual cubic lattice) the interesting pseudo-critical 
region is at very low temperatures. In this region sensible numerical sim- 
ulations are extremely demanding in computer time, due to the extreme 
difficulty in crossing even small barriers. We also believe that a systematic 
comparison of results obtained with different Hamiltonians may be useful in 
finding out those universal features that are independent from the detailed 
form of the Hamiltonian. 

In section (2) we define the model we use and the quantities we measure. 
In section (3) we present the results obtained by using finite size scaling on 
small lattices (from 4 3 to 14 3 ) while in section (4) we present the results 
obtained on a large lattice 64 2 x 128. Finally in section (5) we present our 
conclusions. 

2 The Model and the Observable Physical 
Quantities 

We consider a three dimensional Ising spin glass model on a body centered 
cubic lattice. In this model the lattice sites are labeled by and integer valued 
three dimensional vector i. The spins are defined on each lattice point and 
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take the values —1 or 1. 

The Hamiltonian of the model (with couplings that can take the three 
values and ±1) is 

H W] = ~^Y; J hk(Ti(7k • (1) 

i,k 

The couplings J may be zero or take randomly a value ±1. In the simplest 
version of the models is different from zero if and only if 

\i-k\ = {{i x -k X f +{iy-kyf + {i z -k z fY < r . (2) 

Different models may be obtained by changing the value of r. In the limit 
r^oowe recover the infinite range SK model, while for r = 1 we define the 
usual short range nearest neighbor model. In this paper we will discuss the 
model with r = 3 1 / 2 , which corresponds to have J ^ when all the following 
three conditions are satisfied: 

\i>x k x \ <J 1, \iy ky\ <J 1; \iz k z \ <J 1 , (3) 

and \i — k\ ^ 0. A crucial parameter in the model is the effective coordination 
number z, which is the number of spins that interact with a given spin (for 
r = 1, z — 6; for r = 3? , z = 26). For large values of z the energy 
is proportional to z 1 ^ 2 . On a Bethe lattice (which is a refined mean field 
approximation) the critical temperature may be computed exactly and one 
finds that 

(z-l)tanh(/W) 2 = l ■ (4) 

In this approximation one finds T^^ e =2.08 and T^~^} = 4.93. One 
difficulty with the original r = 1 model is that the hypothetical critical 
temperature is small (about 1.1). Since under a single spin-flip the minimum 
change of the energy is 4, such a low value of the critical temperature implies 
a very small acceptance rate (about 2%) for Monte Carlo steps in which we 
try to change the energy. This effect should disappear for the r = 3 1//2 theory. 
Moreover a different form of the lattice action may be useful to disentangle 
the lattice artifacts from the universal behavior. 

In the particular case of the 3d Ising spin glass a large value of z should 
increase the system reminiscence of to the infinite range model. In a system at 
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the lower critical dimension for high values of z we should see a sharp change 
of behavior from the predictions of the mean field theory to the asymptotic 
low energy behavior. 

In order to define interesting observable quantities it is convenient to 
consider two replicas of the same system (a and r). The total Hamiltonian 
reads 



H = H[a] + H[r] . (5) 
For the two replica system we can define the overlap 

Qi = ViTi , (6) 

which will play a crucial role in our analysis. We will introduce the correlation 
function of two g's as follows: 

G (i) = ^{Qi+kQk) = J2( a i+k a k) (Ti+kTk) = ^2(cr i+ kcrk) 2 • (7) 

k k k 

We can use this correlation function to define a correlation length. From 
high temperature diagram analysis (or from the field theoretical approach) 
we expect that for large separation 

G(i) ~ e -^- . (8) 
We can define an effective mass as 

We expect that at large i 

C 1 = lim m(i) . (10) 

In our numerical simulations we have not measured the full G(i). We 
have measured the zero 2-momentum Green functions 

G^(d) = -±- £ G(i) (11) 

x y y plane 
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where d runs now only in one lattice direction. We will label with a su- 
perscript (°) this kind of quantities. We have also measured the site-site 
correlation function, but only summing over contributions where one single 
coordinate change (by swapping the lattice in a single chosen direction). Here 
the coordinate increment has the form (x, 0,0). We will denote quantities 
defined in this way with a superscript ^\ 

In a similar way in a finite volume we can introduce the quantity 

9=iEft- ( 12 ) 

i 

In the infinite volume limit the spin-glass susceptibility is defined as 



Xo = Jim V < q 2 > , (13) 

where the upper bar denotes the average over the different choices of the 
disorder. 

We expect the spin glass susceptibility and the correlation length to di- 
verge at the critical temperature with the critical exponent 7 and v respec- 
tively. Below the critical temperature in the mean field approach Xo is pro- 
portional to the volume. More generally in the broken replica approach one 
finds that 

lim V <q m >= f 1 dx q(x) m , (14) 

V^-co Jo 

where q(x) is the order parameter function defined in ref. [1]. 

In the high temperature phase no interesting physical predictions can be 
obtained for the usual magnetic susceptibility (divided by (3) defined as 



X = lim V < m? > , (15) 

m being the total instantaneous magnetization (m = yj^i^i)- Gauge in- 
variance implies that at thermal equilibrium 

X=l. (16) 

At T < T c this equality is valid after summing over all configurations with 
the correct Boltzmann weight. If we restrict the sum only to configurations 
in a given equilibrium state this identity does not apply. 
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3 Finite Size Scaling 



We will discuss here results obtained on small lattice sizes, in situations where 
typically L » £. Since our goal is to establish or disprove the existence of 
a critical behavior for T > let us start by sketching the predictions of a 
finite-size scaling analysis. If scaling is satisfied in the vicinity of a critical 
point (at T > 0), we expect 



where 77 are the anomalous dimensions of the operator q defined in (12) and 
£ is the correlation length that is expected to diverge at the critical temper- 
ature. Moreover, to establish the existence of a finite critical temperature 
it is useful to use the Binder parameter to locate the transition point. It is 
defined by 



If a finite T phase transition exists we expect the curves g(T) obtained 
for different lattice sizes to cross (asymptotically for large enough lattices) 
at T c . This is quite a precise method to find the location of a critical point. 
For a T = singularity the same curves will merge in a single curve as 
T — > + . We will see that the possibility that the exponent v characterizing 
the divergence of the correlation length is greater than 1 makes arduous to 
distinguish between these two cases. 

As we have already discussed we want to distinguish among two different 
scenarios. In one case there is a finite temperature transition and the corre- 
lation length diverges like £ ~ (T — T c )~ v . In our finite size scaling analysis 
we will use the large lattice best fit to T c , 7 and v from section (4). If a 
transition exists we have a precise determination of the critical exponents 
and parameters. 

We should note here that if three is the lower critical dimension and 
we have a T = singularity, it is not clear that the scaling relation (17) 
is satisfied. As we will discuss our results suggest that if the scenario of a 
T = phase transition holds such scaling behavior could not hold. This 
violation of scaling appears in the Heisenberg model in two dimensions and 



x ~£ 2 ~V(t), 



(17) 




«<7 2 » 




(18) 
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is a consequence of the existence of the Goldstone modes. In the O(N) 
symmetric Heisenberg model for N > 2 the correct scaling laws contains an 
effective exponent: 

Xo - L>~^ /(|) , (19) 

where 77(0) = 0. The dependence of the exponent on j is due to the instability 
of the T = fixed point. In the N = 2 case, there is no renormalization of the 
coupling constant (i.e. of the temperature). In the low temperature phase 
one gets the simpler equation 

Xo - /(|) , (20) 

where the function rj(T) is not an universal function. Its value at the tran- 
sition point, i.e., r)(T c ), is universal and it is equal to \. 

We have simulated lattices with linear size L = 4, 6, 8, 10, 12 from T = 5.4 
down to the lowest temperature in which we were sure to have thermalized 
(T = 2.6 for L = 4 and T = 3.6 for L = 12). We have computed the overlap 
among two identical copies of the system, defined in (12). 

We have been careful in checking that we have really reached thermal 
equilibrium. We have used as a basic criterion to check that (q) was com- 
patible with zero for each sample. 

We show in fig. (1) the Binder parameter defined in eq. (18), for different 
values of T. We cannot distinguish any crossing, but we better see some 
merging of the different curves. 

From the large lattice results (see section (4)) we can use the values 
T c = 3.27 and 7 = 2.4 (see (36)) and rj ~ (see (36)) to check the consistency 
of the finite size behavior with a finite T transition. In fig. (2) we plot j% 

versus (T — T c ) . The data collapse on a single curve, showing a good 
scaling behavior, on both sides of T c . It is already clear from these first 
data (illustrated in figures (1) and (2)) that it will be exceedingly difficult 
to distinguish between the two candidate critical (with T c = or T c 7^ 0) 
behaviors. 

To understand better what is happening in the pseudo-critical region, for 
T close to 3.3, it is interesting to apply to the model a magnetic field h. We 

2 

expect q to scale as h*. 8 is related to 77 by the hyper-scaling relation 
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Figure 1: The Binder parameter g(T) as a function of the temperature T for 
different lattice sizes. 
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Figure 2: The scaled overlap susceptibility *f versus the scaled reduced 
critical temperature (T — T c ) L». 
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$ = ¥4^- (2i) 

In presence of h the correct definition of the overlap susceptibility requires 
subtraction of the connected part, i.e., 



Xo = Jim V((q*) -(g)). (22) 

For a finite T phase transition the scaling relation (17) is still satisfied, but 
now (we are sitting at T c ) £ diverges like 

£(h) ~ /T 2 ^ . (23) 

(23) only depends on the critical exponent rj. Once we have measured T c , 
and established that a finite T phase transition exists, we can use (23) to 
find r\. 

It turns out that the correct overlap susceptibility we have just defined 
in (22) is not a good observable for checking scaling. It depends on the 
first moment of (q) that is affected by strong finite size corrections. This is 
because the region of negative overlaps with q < is only suppressed in the 
infinite size limit. We have found preferable to study the behavior of the 
non-subtracted (g 2 ), i.e., of the overlap susceptibility defined in absence of h 
divided times the volume. Here we expect the scaling (17) divided times L 3 , 
i.e., a scaling with L with the power —(1 + rj). 

We have run numerical simulations in presence of a magnetic field. In 
figure (3) we show (q 2 ) for several lattice sizes L = 5, 7, 9, 11, 13 and different 
values of the magnetic field (ranging from h — up to h — 1.5). Again 
we find consistency with T c = 3.27. The preferred value for rj is negative 
and close to —0.1. Let us stress that all the finite size scaling fits are not 
giving very precise predictions. There are many free parameters, and that 
makes the fitting procedure questionable. Still we should note that all the 
exponents we find, when assuming a finite T transition, are fully compatible 
with the ones found for the r = 1 model in the previous work of references 
[2, 3, 4, 5]. 

As we already hinted the finite size scaling results are also compatible 
with a T = singularity. We will use the best value (45) of the parameters 
defined in (44). In figure (4) we show the rescaled susceptibility Xo (again 
without magnetic field, now) for the different lattice sizes. The curves for 
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Figure 3: The rescaled (q 2 ) as a function of the rescaled applied magnetic 
field, for different lattice sizes. 
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Figure 4: The rescaled Xo as a function of j. This scaling is compatible with 
a singularity at T = 0. 

different lattice sizes scale tremendously well, and the comparison with figure 
(2) is instructive. This is, as we will discuss in the following, fully compatible 
with the results obtained for the large lattice size, in a regime where £ >> L. 

If the transition is at T = the usual scaling laws imply that the corre- 
lation function at large distance behaves as x~^, with ( = d — 2 + rj. When 
the ground state is not degenerate the T = correlation function goes to 
a constant value at large distance, implying ( = and in 3d r] — — 1. The 
value we estimate for rj turns out to be not so close to —1, and using rj = — 1 
does not make our curves to scale. 

Here we see two options. One possibility is that ( ^ in 3c? Ising spin 
glasses (our best fit is close to ( ~ .6). This possibility cannot be excluded. 
For example in 2d [5] ( is estimated to be in the range .2 ~ .3. In our 
case, where the coupling constants J take the values ±1, the ground state is 
highly degenerate, and there are no general a priori reasons for ( = to hold 
(however it has been suggested in [13] that at the lower critical dimension we 
expect indeed ( = 0). The other possibility is that to get good scaling for rj 
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Figure 5: The overlap susceptibility divided times the asymptotic large lattice 
value (which we denote here by xl) versus j. The line is the best fit to the 
form (31). 

we have to go at lower values of T. Here we have been obliged to seat at not 
too low T's and it is quite possible that the value of eta in this temperature 
range is quite different from its zero temperature limit. 

In figure (5) we have tried to show the scaling behavior in a suggestive 
form. We plot Xo j- L \ as a function of LIE for the different lattice sizes. The 
values of Xo(oo) and £ are those discussed in the next sections and computed 
on very large lattices (which we judge to be free from systematic errors in 
our statistical precision). The data smoothly collapse on a single curve. 

From these data it is not clear if the 3d Ising spin glass undergoes a finite T 
phase transition (and mainly the puzzling behavior of Binder cumulant seems 
to point toward something different). If we assume a finite T c our predictions 
for the critical exponents agree with those reported in the literature (for the 
first neighbor cubic lattice model). 

Though high temperature expansions predict a finite temperature transi- 
tion (which agrees with that found in numerical simulations) we consider the 
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compatibility of our data with a T = phase transition serious (and we will 
discuss this kind of evidence in more detail in next section, when discussing 
our large lattice results). 

As we have already remarked the behavior of the Binder cumulant below 
T c is different from what happens in normal spin systems. It is also very 
different from what one measures in spin glasses in high dimensions, and a 
few more comments are maybe in order. Let us consider what happens in 
the usual ferromagnetic Ising case, by defining the function 

9{T) = hm g(T, L) , (24) 

where here g(T, L) is defined in terms of the moments of the order parameter 
m, the total magnetization of the system. In this non-disordered case we 
have that 



for T <T C g(T) = 1 , for T > T c g(T) = . (25) 

Moreover the quantity g c = g(T c ) is a function of the dimensionality of the 
system. It increases when the dimension decreases, and goes to 1 at the lower 
critical dimension. 

The situation is different in spin glass models. In this case in the mean 
field approximation g(T) is not trivial at low temperature. One finds that 
below T c 

= 3 _ 1 Jdxqjx)* = 3 _ 1 JdP(q)q 4 
2 2{jdxq{xff 2 2(fdP(q)q*) 2 ' 

Using the mean field expression for the dependence of q(x) over T one 
finds that 



g_ = hm g(T) = 1 , (27) 

but the function g(T) is non trivial. The statement g ^ 1 coincides with the 
fact that the P(q) is not equal to a 5-function, and implies replica symmetry 
breaking. In the mean field approximation no closed formula exists for g, 
however one finds that qualitatively g behaves as 

T T 

1-A-(1--). (28) 
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In other words g(T) vanishes linearly both at zero temperature and at the 
critical temperature. For T > T c one still finds that g(T) = 0. Below the 
upper critical dimension (d^ = 6) according to the prediction of ref. [13] 
(?_ becomes different from 1. Slightly below d = 6 the function g(T) is not 
monotonous, but it is possible that it becomes monotonous at sufficient low 
dimensions, i.e., near three dimensions. It is tempting to conjecture that near 
the critical dimension one finds that g c becomes close to g_. It is difficult to 
assess quantitatively the values of these two quantities. If we use our best 
estimate for T c we find g- ~ (.65 ± .05), and a very similar value for g c . We 
can only tentatively conclude that: 

• The L independence of g(L,T) in the (pseudo)-low temperature phase 
and the fact that g(L, T) is different from 1 is a clear signal that replica 
symmetry is effectively broken in this region. Obviously if there is no 
finite T phase transition this symmetry breaking will eventually disap- 
pear for very large lattices, but it will correctly describe the physics 
of the system for large lattices with L smaller than the exponentially 
large correlation length £. 

• The shape of the function g(T) is in qualitative agreement with the 
predictions of the renormalization group and it suggests that the lower 
critical dimension is close to 3 (and very probably exactly 3). 

Let us now discuss in some detail the form of finite size effects. This is 
very interesting, mainly since we have to plan larger scale numerical simula- 
tions, and we want to be sure to optimize the use of our computer time. We 
will describe here the strategy that should eventually lead us to a numerical 
simulation in which we can establish in a clear way which kind of singularity 
the 3d Ising spin glass undergoes. For lattice sizes much larger than the cor- 
relation length one finds that (in presence of periodic boundary conditions) 
the finite volume corrections are exponentially small. The leading correction 
can be computed in perturbation theory, giving: 

Xo(L)= X o(oo)(l-Ce^e^) , (29) 

where C is some computable constant, and A is the coupling constant of a 
3 -like interaction in a field theoretical framework. Close to the critical point 
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the usual scaling laws imply that quantity £ 3 A 2 goes to a constant. So we 
obtain: 



Xo (L) = Xo(oo)(l - C e~ + 0(e— )) . (30) 

We have fitted our data for the correlation length on small lattices, divided 
over the large lattice result, as 

The best fit works very well. We show it in figure (5). For a finite 
temperature transition C is important. It is universal and in principle it can 
be computed in a field theoretical renormalization approach. 

These data are relevant since they are crucial for planning simulations 
free of finite size effects on large lattices. We see that if we require finite size 
effects to be smaller than 1% we need to have L/£ > 6, while to reach a 10% 
accuracy we can accept L/£ > 3.5. 

In a similar way it is interesting to compute 

K = x °_*° . (32) 

Xo 

The quantity K measures the susceptibility system to system fluctuations. 
We expect it to have similar properties to the Binder cumulant g. In partic- 
ular at low temperatures mean field predicts that 

1/ fdxqjx)' x 
K ~ 3\(5 dx q(x)*)* V- { ] 

In other words mean field theory predicts that 



<g 4 > - (g 2 ) 2 = 2((g 2 ) 2 - (q^) ■ (34) 

The size dependence of K can be used to estimate the number of different 
realizations of the quenched disorder we need to extract an accurate value of 

Xo- 

The measurament of K is rather delicate because for each system we 
must know the value of (q 2 ) with high accuracy. In figure (6) we plot K as 
a function of h for L = 6. 
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Figure 6: K as a function of j. 

4 Large Lattice Results and Discussion 

Our large lattice runs have been done on a 64 x 64 x 128 lattice, on the 8192 
processor DECmpp at Syracuse NPAC. We have always studied the evolution 
of two replica of the system in the same realization of the quenched disorder. 
In this way we have been able to compute the overlap between two replica. 

We have studied the behavior of the system for two different realizations 
of the quenched random couplings. We give in table (1) the details about 
the two series of runs (the number of millions of sweeps is for each of the two 
replica we studied in a given coupling realization). 

We studied two different realizations of the random noise mainly to check 
the size of the fluctuations of Xo- We wanted to be sure that even for our 
T point closer to criticality (T = 3.6) sample to sample fluctuations are not 
too dramatic. In fig. (7) we show that in the worst case the two results for 
Xo deviate of less then two standard deviation (in this and in the following 
figures the smooth lines just join the Monte Carlo data points with straight 
segments). But we know from our binning analysis that the error we quote 
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T 


Sample 1 


Sample 2 


6.0 -> 4.4 


0.005 + 0.5 




4.3 -> 3.8 


0.5 + 2.0 


0.5 + 2.0 


3.7 


0.5 + 4.5 


0.5 + 14.5 


3.6 


2.5 + 9.0 


2.5 + 30.0 



Table 1: For the two different realizations of the couplings, the number of 
millions of MC sweeps we used. We give the number of thermalization sweeps, 
plus (+) the number of sweeps used for measuring. 

is probably slightly underestimated at the lower T values. So we find this 
result reassuring, consistent with the serious critical slowing down that we 
are observing and with critical fluctuations. 

The internal energies of the two systems are completely compatible (fig. 
(8)), as it is the specific heat (which we measure both from equilibrium 
fluctuations and from the T derivative of the internal energy, fig. (9)). We 
feel confident that on the 64 x 64 x 128 lattice results do not vary much with 
the sample, and in the following we will discuss results averaged over the two 
realizations of the quenched disorder. 

We have estimated statistical errors by a binning analysis. We have sys- 
tematically blocked the data in coarser and coarser sub-samples, to check 
statistical independence of the configuration groups eventually used for the 
final error analysis. Always but for the two lower T values (3.6 and 3.7) we 
have reached a very reliable estimate of the true statistical error. In the two 
last cases the error seems stabilizing under binning, but the evidence is less 
compelling, and we would allow for a possible small underestimation of the 
statistical error (of less, say, than 50%). 

For T going from 6.0 down to 4.4 we present errors based on 9 blocks of 
order of 50, 000 configurations (the actual measurements were taken just once 
in 200 sweeps). From 4.3 down to 3.8 we have 9 blocks of order 400, 000 con- 
figurations each. At T = 3.7 we have used 5 blocks of 3 x 10 6 configurations, 
and at T = 3.6 6 groups of 6.5 x 10 6 configurations. 

In fig. (10) we plot the final overlap susceptibility, averaged over the two 
coupling realizations, as a function of the temperature T. 

Our main goal has been trying to establish (or disprove) the existence of 
a finite T phase transition for the 3d spin glass model under study. Since 
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0.84 





3.5 4 4.5 5 5.5 6 

T 

Figure 9: As in fig. (7), but the specific heats. The point with a large sta- 
tistical noise are from the energy fluctuations, while the ones with a smaller 
noise are from T derivatives of the internal energy. 
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correlation times diverge very fast when approaching the low temperature 
region (or T c + , if it exists), we are not in an easy situation. On a large lattice 
we have to look at data far away in the warm phase (the one we can check 
and trust have thermalized) , and try to decide which kind of critical behavior 
they have. 

At first we have tried fitting \o with a power divergence at the critical 
temperature T c , i.e., 

A x 

X ^ 1 + W^' (35) 

where the subscript p stands for power fit. We show in fig. (11a) our best fit, 
obtained by using all the data points shown in the figure. The results are 

A* = 19.3 ± 1.1, T c = 3.27 ± .02, 7 = 2.43 ± .05 . (36) 

We do not attach much significance to the statistical errors quoted here. They 
are reasonable estimates of a standard fitting routine, but not the result of 
a detailed study of a very complex 3-parameter fit. We will see in a moment 
that the main issue here is not the statistical error, but the systematic error, 
which is, as far as we can judge from the present data, infinite (see later). 

Obviously one would like to select a T region that would allow exposing 
a good scaling behavior (and to be obliged from the fit to discard a high T 
region where scaling corrections are important and a region close to T c where 
finite size effects become sizeable). This would amount, in some sense, to find 
at least the size of the first corrections to scaling. In the present case we have 
to compromise on the quality of the results in (36), which is, still, reasonably 
good. We have checked that by fitting only points close to T c we get results 
that are not so different from the ones given in (36). For example if we fit 
from T = 5.0 down to T = 3.6 we obtain 7 = 2.67 ± .06, and T c = 3.20. 

Let us repeat that here the problem will turn to be mainly the systematic 
error. 

The second functional behavior we have tried assumes no critical point, 
but an essential singularity at T = 0. We have first tried the form 

Xo ~A*(e ( ^ P -l) + C* , (37) 

where the subscript e stands for exponential fit. The power P turned out to 
be very close to 4 (also for the exponential fit to the correlation length £, see 
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later). We have tried fits with different fixed power P, and for the fit to Xo 
(the fit to has a larger indetermination, see later) we find that a power of 
3 or 5 gives clear worse results than a power 4. So we have eventually used 
the 3-parameter fit to the form 

Xo * A*(e^ 4 - 1) + C* , (38) 

which gives results 

A x e = 1.67 ± .05, Bl = 5.38 ± .01 , C e x = 1.28 ± .05, (39) 

The best fit is very good, and we show it in fig. (lib), on the right. The x 2 is 
much better than for the power fit (12 versus 29 with some slightly arbitrary 
normalization). 

The divergence of the correlation length as a function of (T — T c ) gives, if 
a phase transition exists, the exponent v. We have repeated here the analysis 
we have discussed for Xo- In fig- (12) we give ^ (we have defined before) as a 
function of T. is always compatible with but has a larger statistical 
error. 

Our estimator for ^ is defined by taking the weighted average of the 
effective mass estimator at distance d 

C w (d) 

*< d > s • < 40 > 

for d going roughly from £ to 2£. In this way we are making systematic 
effects (coming from small distance contributions) and statistical error small. 
A typical fitting window is d from 2 to 3 at large T down for example to 8 
to 15 at T = 3.7. We have estimated errors by using a standard binning 
plus jack-knife procedure. Our conclusions about the statistical significance 
of the sample coincide with the ones we have drawn for Xo- 

Also in this case we have tried a power fit and an exponential fit. For the 
power fit we used the form 

(T - T c ) 

with the result 

Al = 2.73 ± .11, T c = 3.24 ± .03, v = 1.20 ± .04 . (42) 
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(41) 




Figure 11: The overlap susceptibility, same than in fig. (10). Here the con- 
tinuous lines are the results of the best fits. In fig. (11a), on the left, the 
result of the power fit to the form (35), and in fig. (lib), on the right, the 
result of the power fit to the form (38^. 
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Figure 12: The correlation length averaged over the two different sam- 
ples. The continuous line is here only to join neighboring points. 
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Even if the results are very reasonable, the fit is not good (as shown in 
fig. (13a), on the left). The x 2 is ver Y high (~ 120), and the points close to 
T c are the one that do not fit (very dangerous caveatl). Still, if we take these 
data seriously, we have to notice that T c is the same we estimated by using 
Xo, and that by means of the scaling relation 

7 = ^(2-7?) (43) 

we get i] ~ 0. 

The exponential fit has the form 



f ~ Ai(e ( ^ 4 - 1) + C\ , (44) 
and gives 

A\ = 1.41 ± .05, Bl = 4.21 ± .02 , C\ = .46 ± .01, . (45) 

Such best fit is very good, and we show it in fig. (13b). The x 2 is 4 times 
smaller than for the power fit. This fit is by far a better fit than the fit to a 
power law behavior. 

For £ the evidence for the power in the exponential being 4 is less com- 
pelling than for Xo- Here fit with power 2, 3 or 5 are acceptable, also if the 
X 2 is minimum at power 4 (or 5, which gives a very similar fit. For power 3 
a small decrease in quality is already apparent). 

In fig. (14) we show the data for 

= Xom 2 , Zi = gim , (46) 

from the data we have already shown for m, the inverse correlation length. 
We expect both quantities to diverge as m v in the small m limit. Both 
quantities can are well fitted with a power law with r\ ~ —.25. 

An independent way to measure 77 is to study directly the data for the 
correlation function G^ 1 ^. At large distances the data can be fitted as 

^e- rar . (47) 

Z(f3) seems to diverge close the critical temperature, with a very small power 
~ m - ' 1 , making this estimate of r\ quite different from the previous one. The 
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Figure 13: The correlation length, same than in fig. (12). Here the continuous 
lines are the results of the best fits. In fig. (13a), on the left, the result of 
the power fit to the form (41), and in fig. (13b), on the right, the result of 
the power fit to the form (44). 
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Figure 14: The quantities Z (lower curve) and Z\ (upper curve) as function 
of m in a logarithmic scale. 
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discrepancy among the two estimates of r] is likely to be related to the small 
asymptotic value of rj. 

As a check we have analyzed the data for the correlation function 




decaying fits to the correlation function are good implies that for s > 1 the 
function C(s) is well approximated by e~ s . At small values of s the function 
should go to zero as s v . Alas, since we cannot reach very small values of s it 
is difficult to use this method to get a precise determination of r\. 

Let us insist on the difficulty in reaching a definite conclusion about the 
critical regime by presenting some more fits (figures (15) and (16). Here we 
are analyzing the overlap susceptibility Xo as function of (3. In fig. (15) we 
show the best fit to the form (35) with the parameters given in (36) (with a 
transition at a critical temperature), and we superimpose a second fit of the 
form 



with A = .085 and B = 15.16. Again, although the two functional forms 
imply a very different critical behavior, in the region we have studied they 
are indistinguishable. 

We can try more. A similar phenomenon is displayed in fig. (16). Here 
we show dependencies that imply a transition at zero temperature: 



In the first best fit we find A = 383 and u = 3.33, while in the second best 
fit we get A = 5.9, B = —69.8 and C = 246. uo turns out to be not so far 
from 4, as we already remarked. 

The four fits all give reasonable results. It is impossible to use the data 
to reject one of them. Of course we could choose the one with smallest x 2 i 
but this procedure may give an incorrect answer since we have neglected 
sub-asymptotic terms, inducing an systematic error which are out of control. 



log(Xo) = A e™ , 



(49) 



log(xo) = Af3 w 

log(x ) = A + BP + CP 2 . 



(50) 
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Figure 15: Two fits to the data for the susceptibility Xo, as a function of (3, 
according to equations (35,49). 
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Figure 16: Two fits to the data for the susceptibility Xo, as a function of (3, 
according to equation (50). 
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From these data, we tend to conclude we have a hint for the absence 
of a phase transition in the 3d spin glass. If on the contrary such a phase 
transition is present, than we have given a reasonably precise estimate of the 
critical exponents. 

5 Conclusions 

We believe we have pointed out an open problem that in recent papers was 
quoted as solved. Nowadays it is usually said that the existence of a phase 
transition is established. For example ref. [7] about aging phenomena (see 
[8] for more aging papers) claims that it is common lore that 3d spin glasses 
undergo a finite T c phase transition. It does not seem to us that the existence 
of a phase transition is well established at all. 

The possibility of 3 being the lower critical dimension is appealing. We 
have in mind a scenario where the predictions of the mean field theory de- 
scribe fairly the behavior of the system down to d — 3, where the transition 
disappears. In no cases, as it is sensible to expect, the system behaves as a 
normal ferromagnet. At low T in 3d the system is reminiscent of the mean 
field picture up to a critical length which is function of T, and diverges at 
T = 0. 

As it was noted many years ago in ref. [14] at the lower critical dimension 
we expect j noise for the power spectrum of the magnetization, that agrees 
with what has been observed experimentally [15]. 

It is clear that there is an apparent critical temperature. Close to this 
pseudo-T c the correlation length becomes so large that it cannot be measured 
on the lattice sizes that are normally studied. Below such temperature the 
system behaves as if it is in the low temperature phase, irrespectively of the 
existence of the transition (think about the Id normal Ising model for low 
values of T). 

The only way to disprove the existence of a transition at finite tempera- 
ture would be to show that the data for the susceptibility and the correlation 
length cannot be fitted with power law singularities at finite temperature. On 
the contrary to present an evidence for a transition at finite temperature one 
should show that the data can be fitted as power law singularities and cannot 
be fitted with functions having only singularities at zero temperature. Our 
data, as well those from the very long simulations of Ogielski and Morgen- 
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Figure 17: T versus r (see the definitions in the text). 



stern [2, 4], can be fitted in both ways. As we already said, we do not think 
that we can discriminate the two admissible behaviors from the value of the 
chi-square x 2 > i- e -> of the quality of the fit, especially in an approach where 
corrections to scaling have been neglected. Unfortunately in absence of clear 
predictions about the low temperature behavior it is difficult to exclude the 
possibility of a transition at T = 0. 

To visually discriminate among the two possibilities we plot in figure (17) 
the quantity 

— dP (51) 



cHog(Xo) 

versus r = log(x ). A finite T transition implies that 



(52) 

with A = I/7, while a transition at T = with a divergence of the form 
exp(/3 w ) implies 
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Figure 18: logXo versus T. Fit 1 (continuous line) is to a power law singu- 
larity, as in eq. (35) fit 2 (dashed line) has a T = singularity, as in eq. 
(49). 



with B — 1 — i. B = 1 corresponds to an behavior. Our best fits 
give A = .29 and B = .86. The data are noisy at high temperature (low 
t). Clearly it is difficult to select one fit, especially since we have neglected 
corrections to scaling. Data seems to prefer a straight line with a coefficient 
not far from one, but we are unwilling to rely on this kind of evidence. 

What can be done with a better numerical simulation? To get a hint 
we have extrapolated two typical fits at a reasonable low T. We show them 
in fig. (18). We have considered a simple power singularity at T ^ 0, and 
a divergence at T = of the form e Ae . From our present best fits we 
can deduce that at, say, T = 3.4, we would be able to discriminate. If the 
data would be really following the finite T singularity scenario (first case), 
the strong increase of the susceptibility could not be fitted by the double 
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exponential scenario, and the zero temperature transition should be refuted. 

In the opposite case, where the hypothetical data would follow a form of 
the second kind (a double exponential singularity at T = 0) we find that a 
power fit would still be a good fit, but with a larger value of 7 and smaller 
value of T c . This variation of the value of the best fit parameters with the 
temperature interval used for the fitting would then be taken as good evidence 
for the existence of a zero temperature transition. 

If the double exponential singularity behavior is correct, the correlation 
length should increase of a factor about 2.5 when going from T = 3.6 to 
T = 3.4. That means that a reliable estimate will be possible on 128 3 lattice, 
only slightly larger of what we used here, and not out of reach of the present 
technology. An increase of the computer time of more than one order of 
magnitude seems unfortunately necessary, but this is also a reasonable goal. 
Such a computation seems possible in a not too far away future. 

It is also possible that a careful analysis of the model at low T could allow 
to show the absence of a phase transition. In this case it would be essential 
to identify the renormalization group flow away from the zero temperature 
fixed point. 
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